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Abstract 

The lepton flavor violation processes J/ip — ► fir and er are searched for using 
a sample of 5.8x10 J/ip events collected with the BESII detector. Zero and one 
candidate events, consistent with the estimated background, are observed in J/ip — ► 
fir, t — ► ev e v T and J/ip — ► er,T — ► iiv^Vr decays, respectively. Upper limits on the 
branching ratios are determined to be Br(J/ip — ► /j,t) < 2.0 x 10~ 6 and Br(J/ip — ► 
er) < 8.3 x 10~ 6 at the 90% confidence level (C.L.). 

PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 



1 Introduction 



In the Standard Model ( SM), lepton flavor is conserved, but it is expected to 
be violated in many extensions of the SM, such as grand unified models [1], 
super symmetric models [2] , left-right symmetric models [3] , and models where 
electroweak symmetry is broken dynamically [4]. Recent experimental results 
from Super-Kamiokande [5], SNO [6], and KamLAND [7] indicate strongly 
that neutrinos have masses and can mix with each other. Consequently, lepton 
flavor symmetry is a broken symmetry. There have been many studies both 
experimentally and theoretically on searching for lepton flavor violating (LFV) 
processes [8], mainly from /i, r and Z decays [9]. Theoretical predictions of 
LFV in decays of charmonium and bottomonium systems are discussed in 
Refs. [10,11,12], and the search for the J/ip — > e/i LFV process at BESII is 
presented in Ref. [13]. In this paper, we report on a search for LFV via the 



decays J/ip — > /ir, r — > ev e v T and J/^ -^ er,r ^ nVpV r using 5.8 x 10 7 J/ip 
events collected with the BESII detector. 



2 BES detector 



The Beijing Spectrometer (BES) [14,15] is a conventional solenoidal magnetic 
detector at the Beijing Electron Positron Collider (BEPC). The upgraded 
version of the BES detector, BESII, includes a 12-layer vertex chamber (VC), 
surrounding the beam pipe and providing trigger information; a fourty-layer 
main drift chamber (MDC), located radially outside the VC and providing 
trajectory and energy loss (dE/dx) information for charged tracks over 85% of 
the total solid angle; and an array of 48 scintillation counters surrounding the 
MDC to measure the time-of-flight (TOF) of charged tracks with a resolution 
of ~ 200 ps for hadrons. The momentum resolution of the MDC is cr p /p = 
1.78%\/1 + p' 2 {p in GeV/c), and the dE/dx resolution for hadron tracks is 
about 8%. Radially outside the TOF system is a 12 radiation length, lead- 
gas barrel shower counter (BSC). This measures the energies of electrons and 
photons with an energy resolution of ue/E = 21%/y/E (E in GeV). Outside 
the solenoidal coil, which provides a 0.4 Tesla magnetic field over the tracking 
volume, is an iron flux return that is instrumented with three double layers of 
counters to identify muons of momentum greater than 0.5 GeV/c. It provides 
coordinate measurements with resolutions in the outermost layer of 10 cm and 
12 cm in r<f) and z. The solid angle coverage of the layers is 67%, 67%, and 
63% of An, respectively. 

In the analysis, a GEANT3 based Monte Carlo program (SIMBES) with de- 
tailed consideration of detector performance (such as dead electronic channels) 
is used. The consistency between data and Monte Carlo has been checked in 
many high purity physics channels, and the agreement is reasonable. 



3 Event selection 



We require candidate events for J/ip —* fir, r — > ev e v T and J/ip — ► er,r — »■ 
liv il y T to have two well reconstructed and oppositely charged tracks, each 
of which is well fitted to a helix originating from the interaction region of 
\x\ <0.015 m, \y\ <0.015 m, and \z\ < 0.15 m and with a polar angle, 9, 
satisfying | cos#| < 0.8. To reject cosmic rays, the time of flight difference of 
the two charged tracks should be less than 4 ns. 

Isolated photons are defined as those photons having an energy deposit in the 
BSC greater than 50 MeV, an angle with any charged track greater than 15°, 



and an angle between the direction defined by the first layer hit in the BSC 
and the developing direction of the cluster in the xy-plane less than 18°. There 
must be no isolated photon in the selected event. 

Information from the BSC, TOF, and MDC (dE/dx) is used to select electrons. 
Fig. 1(a) shows the ratio of the energy deposited by the electron in the BSC to 
its momentum (E/P) for Monte Carlo simulated events, and Fig. 1(b) shows 
the energy deposited by the muon in the BSC for Monte-Carlo simulated 
events. To be an electron, the charged track should have no hits in the muon 
counter, and the E/P ratio should be larger than 0.7. To further distinguish 
the electron from hadrons, it is required that p e dE/dx > p n dE/dx , p e dE/dx > p^ E/dx 
and Pxof > Ptof? where p dE / dx and p l TOF are the particle identification 
confidence levels for the dE/dx and TOF measurements and % denotes e, ir, K 
or p. 
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Fig. 1. (a.) Distribution of E/P for electrons (MC simulation), (b.) Distribution of 
energy deposited by muons in the BSC (MC simulation). The solid histogram repre- 
sents J/tp — ► er,r — > iiv^Vt channel, and the dashed one is for J/ip — > /xr, r — > ev e v T 
channel. 



To select muon tracks, the differences, 5i(i = r<f), z), between the closest muon 
hit and the projected MDC track in each layer are used. A good hit in the \i 
counter requires \Si\ < 2<7j for i — r<f> and z. The total number of good \x hits 
in the jx counter, nf°° , can range from to 3. A track is considered as a muon 
if the deposited energy in the BSC, shown in Fig. 1(b), is less than 0.3 GeV 
and /j, 9 f°° is equal to 3. 

For the decay of J/ip — *■ er,r — >■ yUZ/ M z/ T , the momentum of the electron is 
monochromatic, while that of the muon is broad, as shown in Fig. 2. The 
main background for this channel comes from J/ip — > ( / ~f)fi + fi~ and e + e~ — > 
(7)/x + /i _ , which is shown as the dashed histogram in Fig. 2. This background 
can be rejected by requiring that the momentum of the electron P e be in the 
region from 1.00 to 1.08 GeV/c and the momentum of the muon be less than 



1.4 GeV/c. Similar requirements P e < 1.4 GeV/c and 1.00 < P M < 1.08 GeV/c 
are applied to J/?/> — > fir, r — ► ez/ e z/ r candidates to suppress the background 
from J /if) — ► (7)e + e~ and e + e~ — > (7)e + e~. Applying these requirements, no 
candidates for J/ip — > fir, r — > eu e u T and one candidate for J ftp — > er,r — > 
liu^u T survive. 
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Fig. 2. Monte-Carlo distributions of muon momentum. The solid histogram repre- 
sents J/tp ^ er,T ^ iw^v-f channel, and the dashed one is for the J/ip — ► (y)(i + [x~ 
channel. 



4 Efficiencies and backgrounds 



In this analysis, the \x particle identification efficiency e^pju in the /i counter is 
determined using real \x tracks. All other efficiencies, including the geometric 
acceptance, momentum requirement efficiency, electron particle identification 
efficiency, etc., are combined into one term, €mc, which is determined by 
Monte-Carlo simulation. The overall efficiency is calculated as ^totai— e ixPiD x 
£mc- 

The \i track sample selected from 5.8 x 10 7 J/ip — ► (•y)fj J + fj,~ decays, as de- 
scribed in Ref. [13], is used to determine the /i particle identification efficiencies 
in both channels. The ji particle identification efficiency is a function of the 
transverse momentum, P xy , of the muon. Therefore, e^pju is determined from 
Y^^i^ii where e$ is the \x particle identification efficiency in the ith P xy bin 

i 

determined from the \x track sample, and Wi is the weight corresponding to the 
number of events in the bin determined from the signal MC. Table 1 lists the e* 
and Wi in the different P xy regions, and Table 2 lists the selection efficiencies. 

The remaining background in both the J/ip — > fir, r — > ev e v T and J/ip — > 
er, r — > [iv^Vt processes are studied through Monte Carlo simulation. Almost 
all two-prong decay modes are generated with 5 to 10 times the number of 



Table 1 

The €i and w% values in different P xy regions in the J ftp — ► jjlt and er channels. 







J/ip — ► //r, r — ► ev e v T 


J/V> ->■ er,r -> /iP^ 


P xy (GeV/c) 


Ci (%) 


Wi 


tJ7j 


0.5 < P^ < 0.6 


0.0 


0.0 


13.7 


0.6 < P xy < 0.7 


0.0 


3.1 


12.3 


0.7 < P xy < 0.8 


8.1 


11.3 


10.9 


0.8 < P xy < 0.9 


40.6 


17.8 


9.7 


0.9 < P xy < 1.0 


52.2 


33.9 


7.9 


1.0 < P X y < 1.1 


53.4 


33.9 


6.3 


1.1 < Psj, < 1.2 


56.2 


0.0 


4.5 


1.2 < Pxj, < 1.3 


57.7 


0.0 


2.8 


1.3 < P Xtf < 1.4 


53.6 


0.0 


1.0 



Table 2 

Efficiency summary 





3/tp — ► /LIT, r — v ei/ e i> T (%) 


J/V» ->■ er, r ->■ /uP^ r (%) 


e^piD 

(-MC 


43.9 
26.2 


17.0 
28.1 


ttotal 


11.5 


4.8 



events expected from 5.8 x 10 7 J/ip events. For J/ip — > er,r — > \iv il y T ^ the 
estimated background is about 0.4 events from <//■?/> — > K* (892)~ K + (+c.c.) . 
For the decay J/ip — ► ^ut, r — ► eP e z/ r , no simulated events survive. 



5 Systematic errors 



The systematic errors in the branching ratio measurements come from the 
uncertainty of the MDC tracking efficiency for charged tracks, the error from 
the number of J/ip events, the differences in the efficiencies between data and 
Monte-Carlo simulation for some selection criteria, such as the electron and 
muon identification criteria, as well as the uncertainty in r decay branching ra- 
tio. The systematic errors from each source are listed in Table 3; the dominant 
error is from muon identification. Adding all the systematic errors in quadra- 
ture, the total systematic errors are 16.9% and 15.4% for J/ip — > fir, r — > ev e v T 
and J/tp — ► er,r — > /iP^Ur respectively. 



Table 3 
Summary of systematic errors 



Source 


J/V- 


-*• flT,T -*• 


ev e 


M%) 


J/Y>- 


-> er 


,T -»■ IJ,9^V T (%) 


epiD 




3.5 










3.3 


HPID 




15.4 










13.7 


Br(r decay) 




0.3 










0.4 


MDC tracking 




4.0 










4.0 


Number of J/ijj events [16] 




4.7 










4.7 


Sum 




16.9 










15.4 



6 Results and discussion 



No J /ip — > /it, r — >■ e^e^T candidate and one J/ip — ► er, r — *■ fiu^Ur candidate 
are observed from a sample of 5.8xl0 7 J/?/> events, where the estimated back- 
ground events in the two chanels are of and 0.4 events, respectively. The 
background events are ignored for the conservative estimation. Upper limits 
on the branching ratios of J/ip — > \xt and J/ip — > er are calculated with: 



Br(J/i) ^X)< 



X(N Signal ,N 



BG 



Nj/^ x Br(X ^Y)x e J/v ,^ x _ 



• Y 



where X and Y stand for the intermediate and final states, A is the upper 
limit on the number of observed events at the 90% C.L., Nsignai and Nbg 
are the numbers of observed signal and background events respectively, Nj/^ 
represents the total number of J/ip events, and e is the detection efficiency. The 
values of \{N Signa i and Nbg) can be calculated using the method described in 
Refs. [17] and [18]. 

With the numbers summarized in Table 4, the upper limits on the branching 
ratios, after incorporating the systematic errors, are 



Br{J/i) - 
Br(J/i[> 



fir) < 2.0 x 10~ 6 , 
> er) < 8.3 x 10~ 6 



at the 90% C.L. 



Previously BES reported an upper limit on Br{J/ip — > e/i) to be 1.1 x 10 6 
at the 90 % C.L. [13]. 

In summary, the LFV processes J/ip — > fir and er are searched for using 
a sample of 5.8 x 10 7 J/ip events. No candidate for J/ip — > fir,r — > ev e v T 



and one candidate for J/ip 



/iz/„z/ T , consistent with the estimated 



Table 4 

Numbers and efficiencies used in the calculation of the upper limits. 





J/V 


-^ flT,T ^ 


ev^r 


J/1> 


->■ er, r ->■ /U^z/ T 


Nbg 














e(%) 




11.5 






4.8 


Nj/i, 




5.8 x 10 7 






5.8 x 10 7 


^Signal 











1 


HNsignal,N bg ) 




2.4 






4.0 


Br(r decay) (%) 




17.84 






17.37 


Upper limit of Br. 




2.0 x 1(T 6 






8.3 x 10~ 6 



background, are observed. The upper limits on the branching ratios at the 
90% C.L. are determined to be Br{J/ip — » /ir) < 2.0 x 10~ 6 and Br(J/ip -> 
er) < 8.3 x 10~ 6 . 
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